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ABSTRACT 

Aims. There are still many open questions as to the physical mechanisms at work in Low Luminosity AGN that accrete in the extreme 
sub-Eddington regime. Simultaneous multi-wavelength studies have been very successful in constraining the properties of Sgr A*, 
the extremely sub-Eddington black hole at the centre of our Milky Way. M81*, the nucleus of the nearby spiral galaxy M81, is an 
ideal source to extend the insights obtained on Sgr A* toward higher luminosity AGN. Here we present observations at 3 and 1 mm 
that were obtained within the framework of a coordinated, multi-wavelength campaign on M81*. 

Methods. The continuum emission from M81* was observed during three epochs with the IRAM Plateau de Bure Interferometer 
simultaneously at wavelengths of 3 and 1 mm. 

Results. We present the first flux measurements of M81* at wavelengths around 1 mm. We find that M81* is a continuously variable 
source with the higher variability observed at the shorter wavelength. Also, the variability at 3 and 1 mm appears to be correlated. 
Like Sgr A*, M81* appears to display the strongest flux density and variability in the mm-to-submm regime. There remains still 
some ambiguity concerning the exact location of the turnover frequency from optically thick to optically thin emission. The observed 
variability time scales point to an upper size limit of the emitting region of the order 25 Schwarzschild radii. 

Conclusions. The data show that M81* is indeed a system with very similar physical properties to Sgr A* and an ideal bridge 
toward high luminosity AGN. The data obtained clearly demonstrate the usefulness and, above all, the necessity of simultaneous 
multi-wavelength observations of LLAGN. 

Key words. LLAGN - mm-observations - galactic nucleus 



?H 1. Introduction which lead to differing predictions on the detailed spectral en- 
^ ergy distribution (SED) and its variability. Another open ques- 
At high accretion rates onto black holes, the infall of gas tion is the role of jets in dissipating accretio n energy (e.g. 
and dissipation of energy via mostly thermal processes m a iBlandford & KonigllTTOTOl: IP^cke & BiermannI fT995h and the 
thin disk are fairly well understood in general terms. In con- extent of their contribution to the SED, particularly at higher fre- 
trast, there is still great debate about which physical mech- gugncies ( e.g. Falcke & Markofl 2000; Markoffetal. 2001b.i 
amsms are dommant m objects which are accreting at well |Yuan et alJl2002l) . Since iets are observed in almost all accreting 
below the Eddington rate. The most extremely sub-Eddington astrophysical systems, the importance of jets in weakly accreting 
source currently accessible to observations that allows still rea- systems is a matter of particular interest 
sonable statistics for modelling the accretion/emission mech- 
anisms is the central supermassive black hole of our Galaxy, with a luminosity -lO"''- -""'xLEdd, Sgr A* is - in terms of 
Sagittarius A* (Sgr A*). Its low luminosity has perplexed theo- Eddington luminosity - the weakest accreting black hole with 
rists for decades, and stimulated the re-emergen ce of radiatively observational statistics good enough to fit models to its spec- 
inefficient accretion fl ow models (RIAFS, see |Meha & Falckg trum. It is the primary testbed for theoretical models of extreme 
|2001|: |Quataeri|2003|), including the advection domin ated ac- sub-Eddington accretion, which rely largely on the available 
cretion flow model (ADAF, e.g. [Narayan & Y||1994|) and its radio /submm/NIR /X-ray observations of Sgr A* (for a review 
derivatives. These models are characterised by low efliciency see iMelia & Falckd 120011) . The recent discovery of quiescent 
in converting thermal energy into electromagnetic radiation, and flari ng emission frorn Sgr A* at both X-ray and NIR wave- 
but differ in terms of tiie inner boun dary conditions such as lengths jBaganoff et alj 120011: iGenzel et all 120031; iGhez et all 



the p resence of strong winds (e.g. |Bl andford & Begdmanl |2004l ) has revolutionised our understanding of this particular 



|1999|) and/or convection dQuataert & Gruzinov 1999), each of Low Luminosity Active Galactic Nucleus (LLAGN). The strin- 

gent constraints from these data have ruled out several models. 

Send offprint requests to: R. Schodel leaving only RIAF models, jet models, and combinations thereof 
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as contenders. Recently, lEckart et alj (|2004 l2006l) were suc- 
cessful in obtaining the first simultaneous measurements of the 
emission from Sgr A* at X-ray /NIR wavelengths (with quasi- 
simultaneous sub-mm observations). These and ongoing coordi- 
nated multi-wavelengths campaigns deliver the decisive obser- 
vations for constraining/eliminating weak accretion models fur- 
ther. 

The nearby spiral galaxy M81 (NGC 3031) is an Sb spiral 
galaxy similar to the MiUcy Way. It is l ocated at a distance of 
3.63 + 0.34 Mpc d Freedman et"al]ll994 . lDevereux et al.1 (l2003h 
used spectroscopic measurements of the Ha + [Nil] emission, 
probably emitted from a rotating gas disk inclined at an angle of 
14° + 2°, to infer a mass of 7.0!j x 10^ M© for the central black 
hole in M81. 

The nucleus of M81, termed M81*, shows typical signs of 
AGN activity. It has a power-law, variable X-ray continuum 
jlshisaki et al.ll996l:lPage et al.l2004l) . The X-ra y flux from M81 
is hig hly variable, at scales from days to years dLa Parola et al.l 
|2004|) . The nucle us of M81 displa ys double peaked, broad Ha 
emission lines (Bower et al. 111996*). However, the overall lumi- 
nosity and AGN characteristics of M81* are rather weak and 
the galaxy is classified as a LINER (lo w-ionisation nuclear 
emission-line region, e.g. iHo et al.1 [19961) . With a luminosity 
of the order of 10^^ erg s ' in the radio and 10 ^° ergs"' in the 
optical/X-ray domains (see, e.g., compilations bv lHo et al.l 19961: 
[Ho 1999), its luminosity is < 10"^ times the Eddington luminos- 
ity in any wavelength regime. M81* is therefore counted among 
the low-luminosity AGN (LLAGN). It shows the typical spectral 
energy distribution (SED) of this class of sources, that is char- 
acterised by the absence of the so-called big blue bump, the ul- 
traviolet excess found commonly in the higher power AGN (e.g. 
lHolll999l) . 

At cm-wavelengths, M8 1 * shows large-amplitude variations 
(factors up to two at 2 cm) with timescales of a few months 
and weaker c hanges of the flux density on timescales < 1 day 
jHo et al.1 [199 9). Multi-epoch VLBI observations of M81* at 
~0.01 pc resolution at 20 epochs over 4.5 yr reveal a stationary 
core with a variable (on timescales of ~ 1 yr) one-side d jet of 
length 1 mas (3600 AU) towards the northeast (Bietenh olz et alj 
[2000 ). As for the polarisation prope rties of M81*, circula r po- 
larisation was_detected_at 4;8j8. 4 (iBrunthaler et alj[200ll) . and 
15 GHz (iBrunthaler et al.l [20061) . while linear polarisation ap- 
pears to be absent at these frequencies. This is an intriguing sim- 
ilarity to Sgr A*, where circular polarisa tion also dominates o ver 
linear polarisation at these wavelengths (iBower et al.lll999allbl[cl : 
lAitken et al.'' 200aiBower eran[2003l[2005l) ! 

[Sakamoto et al.l (1200 Ih present observations of the central 
kiloparsec of M8 1 at a wavelength of 3 mm in the CO 7 = 1-0 
line and continuum at lOOpc resolution. They detect molecular 
gas in a pseudoring or spiral arm at about 500 pc, but no giant 
molecular cloud within about 300 pc of the nucleus. They find 
significant intraday variation of the continuum emission from 
M8 1*, su ggesting an emitting region of ~100 AU. 

iReuter & Lesch C1996.) obtained a spectrum of M8 1 * from 
the radio to the mm-regime. They find an inverted spectrum 
up to 100 GHz. Its flux density can be described well by the 



law 5v 



,1/3 



exp -(v/vc) with a turnover frequency of 



200 GHz. They point out the remarkable similarity between the 
spectrum of M8 1 * and Sgr A* and conclude that the same physi- 
cal mechanisms might operate in both galactic nuclei. Therefore 
iReuter & LeschI d 19961) suggest that M8 1 * may be a by ~ 1 0"* up- 
scaled version of Sgr A*. 

Thus M8 1 * is a unique source for comparison with Sgr A* 
and more powerful AGN, and constitutes the next logical step 



after the successful multi- wavelength observations of Sgr A*. 
With this aim, a coordinated, multi-wavelength campaign took 
place in the first half of 2005, involving instruments from the X- 
ray to the radio domain: the Chandra X-ray observatory, the Lick 
telescope (NIR), the SMA, the Plateau de Bure Interferometer 
(PdBI), the VLA, and the GMRT. A compilation of all observa- 
tions and the interpretation of the multi-wavelength data is pre- 
sented by Markoff et al. (in preparation). A detailed description 
of all instruments and the related data reduction would overload 
the multi-wavelength paper. Therefore, some papers are dedi- 
cated to the observations with specific instruments. The X-ray 
emission lines measured in the Chandra observations are dis- 
cussed by Young et al. (in preparation). In this paper we focus 
on the three epochs of mm-observations of M8 1 * that were ob- 
tained with the PdBI during this campaign. 

2. Observations and data reduction 

M81* was observed with the PdBI on 24 February, 14-15 July, 
and 19-20 July 2005. Six antennae were used for the February 
observations in the B configuration of the PdBI, which provides 
typical beam sizes around 1.5" at ~100GHz and around 0.8" at 
~200GHz. The July observations were done with five antennae 
in the more compact D configuration, with typical beam sizes 
around 5" at -lOOGHz and 2.5" at ~200GHz. AU observations 
are listed in Table[T] 

The observations aimed to detect continuum emission from 
M8 1 *. However, due to good conditions the receivers were tuned 
to the '^CO 7=1-0 and 7 = 2-1 transitions at 115.3 and 
230.5 GHz in the February observations. Thus, it was possible 
to search for compact CO-emission within -20" of the nucleus, 
while the continuum could be extracted from the line-free chan- 
nels (The results of the CO-line imaging will be discussed in a 
forthcoming paper). For greater phase stability in the July obser- 
vations, an epoch during which the atmosphere at the site con- 
tains more water vapour, the receivers were tuned to frequencies 
of 80.5 and 86.2 GHz at 3 mm and 24 1 .4 and 2 1 8 .2 GHz at 1 mm. 

The sources 3C273 (February), 3C454.3 (July 14-15) and 
1741-038 (July 19-20) were used for bandpass calibration. 
Phase calibration was performed with the sources 1044+719 and 
0836+710. Primary flux calibrators to determine the efficiencies 
of the antennae were the sources 1044+719 for February 23- 
24 (1.6Jy at 3mm/l.lJy at 1mm), 1044+719 (1.8 Jy at 3mm) 
and 2200+420 (8.7 Jy at 1 mm) for July 14/15, and MWC349 
(1.0 Jy at 3 mm) and 3C454.3 (33.0 Jy at 1 mm) for July 20. The 
phase calibrators 1044+719 and 0836+710 were used to fit the 
time-dependent fluctuations of the amplitude for all baselines. 
Various tests were performed for estimating the uncertainty of 
the flux calibration, e.g. using different primary flux calibrators 
and comparing the resulting fluxes of all observed sources. As an 
additional test, from a comparison of the calibrated fluxes of the 
sources 0836+710, 1044+719, and MWC349 between the three 
epochs we estimated the uncertainty of the absolute flux calibra- 
tion. To provide a quantitative measurement of the data quality. 
Table [2] lists the ranges of the rms values of the phase and am- 
plitude that were obtained during calibration of these values for 
the different baselines for each observing epoch. The table also 
provides values for the estimated uncertainty of the absolute flux 
calibration at 3 mm and 1 mm. As can be seen, the quality of the 
February data is highest. The data from 20 July are clearly better 
than the data from 14/15 July. 

Individual scans of 20 min duration were extracted from the 
calibrated data. Subsequently, a point source was fitted to the re- 
sulting UV tables in order to determine the flux of M8 1 * and of 
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Table 1. Observations of M81* with the PdBI during 2005. V3mm is the exact frequency used around 3 mm, vimm the one at a 
wavelength of 1 mm. Nam is the number of antennae, "Config" refers to the antenna configuration used. 



Start [UT] 


End [UT] 


V3mm [GHz] 


Vl,nn, [GHz] 


NAnt 


Config 


24-FEB-2005 01:11 


24-FEB-2005 19:45 


115.3 


230.5 


6 


Bp 


14-JUL-2005 06:50 


15-JUL-2005 13:51 


80.5 


241.4 


5 


D 


19-JUL-2005 23:17 


20-JUL-2005 16:07 


86.2 


218.2 


5 


D 



Table 2. Data quality: The table lists the ranges of the rms values of the phase and amplitude between all baselines for each observing 
epoch and wavelength. Also, the estimated systematic error of the absolute flux calibration is indicated. 





3 mm 


3 mm 


3 mm 


1 mm 


1 mm 


1 mm 


Date 


phase rms [deg] 


amp rms [%] 


absolute 


phase rms [deg] 


amp rms [%] 


absolute 


24-FEB-2005 


10-20 


5-8 


<10% 


20-40 


11-18 


<15% 


14/15-JUL-2005 


15-35 


8-11 


<15% 


40-60 


25-30 


<30% 


19/20-JUL-2005 


12-28 


~3 


<15% 


30-70 


11-18 


<20% 



the two calibrators 0836-1-710 and 1044+719. Fitting Gaussian 
functions to the UV tables of MS 1 gave very similar results, with 
deviations < 5% at 3 mm and < 15% at 1 mm. Due to the larger 
number of parameters the uncertainties are larger when fitting 
Gaussians. The general shape of the radio light curves is not af- 
fected by the model that is chosen to fit the data. Since MS 1 * 
is known to be a point-like radio continuum s ource that shows 
struct ure only on scales of 1 milli-arcsecond ( iBietenholz et al.l 
I2000 h. we decided that fitting a point source appears to be the 
appropriate way to measure the flux density of MSI*. 

3. Variability of M8V at mm wavelengths 

Figure [U shows the observed millimetre light curves of MSI* 
along with the corresponding data for the two phase calibrators. 
The scaling has been chosen such that the average flux of each 
source equals a value of 1 (arbitrary units), with the fluxes of 
the calibrators shifted for better comparison. Only the error bars 
of the relative uncertainty are shown. As for the absolute uncer- 
tainty of calibration, it would merely shift the curves along the 
vertical axes and have no influence on the shape of the curves 
under stable observing conditions. The flux of the calibrator 
sources should then fluctuate randomly and with small standard 
deviation around a constant value. This is the case for the 3 mm- 
data on 24 February and 20 July, while the conditions were less 
favourable on 14/15 July. In case of the 1 mm-data, the calibra- 
tor measurements fluctuate around a stable average on 20 July, 
they show trends in parts of the data from 24 February (that are, 
however, significantly smaller and significantly less pronounced 
than the trend seen in the flux from MSI*), and show significant 
deviations on 14/15 July. 

As for the overall quality of the data, it can be best asserted 
when examining the calibrator sources (see also Tab.|2]for more 
quantitative information on the data quality). Clearly, the qual- 
ity of the data is higher at 3 mm than at 1 mm due to the re- 
duced phase stability at the shorter wavelength. The best data 
were obtained on February 24, when the lower winter tempera- 
tures and vapour content of the atmosphere facilitated measure- 
ments at 1 mm. The July 20 observations are excellent at 3 mm 
and still good at 1 mm. Observing conditions were highly vari- 
able on 14/15 July, but it was still possible to obtain good data at 
3 mm and acceptable ones at 1 mm. However, a large part of the 
data had to be discarded during calibration so that in total only 
about 6 hours of usable data were left. 

Some variability of the calibrators can be seen in Fig.[T]that 
may be correlated with the variability of the target. We tried to 



take this possibility into account by introducing an additional 
calibration step that forces the flux of the calibrators to be con- 
stant (on average). Thus, any possible systematic uncertainties, 
leading to spurious variability, can be removed. Figures |2] to |4] 
and all values in Tab.[3]refer to measurements after this final cal- 
ibration step was applied (uncertainties introduced by the cali- 
bration process were taken into account). Generally, the effect 
of this additional calibration step is minor It is illustrated for the 
Feb. 24 data in Fig.|2] Since (quasi)-simultaneous measurements 
of sufficient quality for both calibrators were not available for all 
measurements of MS 1 *, there are less data points in some of the 
plots in Fig. [3] and ID than in Fig.[T] 

The calibrated flux density measurements of MS 1 * at 3 mm 
and 1 mm are shown in Fig. [3] (only error bars of the relative 
uncertainty are shown) and the corresponding spectral indices 
in Fig. ID (here, the absolute uncertainties of the flux calibration 
have been taken into account additionally). The spectral index 
is defined by 5v oc v", where Sy is the flux density at a given 
frequency v. Time averages of the flux densities and spectral in- 
dices are listed in Table [3] 

We find the following characteristics of the variability of 
MS 1 * at wavelengths of 3 mm and 1 mm: 

- Variability at the two wavelengths appears to be correlated. 

- The amplitude of variability is by a factor ~1.5 larger at 
1 mm. 

- On 14 July, the flux density at 3 mm is a factor of ~ 3, at 1 mm 
a factor of ~2 higher than at the other two epochs. 

- The spectral index varies between ~0 and — 0.5 (Fig. |4] and 
Tab,[21). 

- MS 1 * shows intraday variability, with the fastest variability 
observed between OS and 12 h UTC on 24 February. 

The light curve of February 24 is shown in the upper panels 
of Fig. [1] The average flux density and its standard deviation is 
indicated for different sections of the light curve. A > 5cr de- 
crease can be seen at 3 mm between 7 and 12 UT. It is accompa- 
nied by a similar decrease at 1 mm. The February 24 light curves 
are shown in Fig. 2 after forcing the two calibrator sources to 
have a constant flux (on average). This additional calibration step 
does not alter the intrinsic shape of the light curves in a signif- 
icant way. Therefore the observed intra-epoch variability must 
have been intrinsic to the source. 

An important question is whether the detected variability 
could be related to polarisation. We can discard this possibility, 
however, because a) The directions of linear polarisation of the 3 
and 1 mm receivers at the antennae of the PdBI are orthogonal. If 
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M81 at 115.3 GHz - 24 Feb 2005 



M81 at 230.5 GHz - 24 Feb 2005 
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M81 at 80.5 GHz - 14/15 Jul 2005 
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MSI at 86.2 GHz - 20 July 2005 
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Fig. 1. Light curves of M81* and of the caHbrators 1044-1-719 and 0836-H710 from the PdBI observations at 3 and 1 mm. The x- 
axes show UT in hours. The fluxes were all scaled to an average value of 1, with the curves of the calibrators shifted for better 
comparison. The vertical lines in the top panels indicate sections for which the average flux and the standard deviation of the 
individual measurements have been calculated. The corresponding values are indicated. 



a sinusoidal pattern due to polarisation were observed at 3 mm, 
one would expect to observe a similar pattern at 1 mm, shifted by 
several hours. We do not detect the expected sinusoidal pattern 
in the light curves nor are the light patterns shifted by several 
hours in time relative to each other, b) Recent observ ations of 
M81* with the BIMA array by lBrunthaler et alJ (l2006l) indicate 



absence of or very low upper limits of a few percent on linear 
and circular polarisation of MS 1 * at 86 and 230 GHz. 
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M81 at 115.3 GHz - 24 Feb 2005 



M81 at 230.5 GHz - 24 Feb 2005 
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Fig. 2. Light curves of M8 1 * and of the calibrators 1044+719 and 0836+710 from the PdBI observations at 3 and 1 mm on February 
24. The light curves were corrected for possible systematic variations by forcing the calibrator sources to have constant flux (on 
average). 



Table 3. Average flux densities of M81*, 5 3 mm and S 1 mm, measured at V3mm and vi mm during the three epochs of the campaign. 
The mean and standard deviation (not the error of the mean) were calculated from the unweighted individual measurements shown 
in Fig. [3] The last column lists the average spectral indices and their standard deviation (not the error of the mean) calculated from 
all individual measurements as shown in the right panels of Fig.|4] 





Stai-t [UT] 


End [UT] 


Vjmm [GHz] 


S 3 mm 


vi ,„„, [GHz] 


S 1 mm 


a 




24-FEB-2005 01:11 
14-JUL-2005 06:50 
19-JUL-2005 23:17 


24-FEB-2005 19:45 
15-JUL-2005 13:51 
20-JUL-2005 16:07 


115.3 
80.5 
86.2 


88.0 ± 11.7 
241.2 ±33.8 
118.7 ± 11.4 


230.5 
241.4 
218.2 


85.6 ± 17.8 
181.2 ±39.1 
74.8 ± 13.3 


-0.06 ±0.16 
-0.38 ±0.17 
-0.51 ±0.15 




M8I at 115.3/230.5 GHz - 24 


Feb 2005 


M81 at 80.5/241.4 GHz - 14/15 Jul 2005 
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1 1 I '11 1 , 



UT [h] 



UT [h] 
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Fig. 3. Flux density of M8 1 * on 24 February, 14/15 July, and 20 July 2005 . The black boxes mark the flux density at 3 mm and the 
green (gray) crosses mark the flux density at 1 mm. The exact frequencies are indicated in the titles of the individual panels. The 
error bars indicate relative, not absolute, uncertainties. 



4. Discussion 

Theory and observat ions indicate that s ub-Eddington black holes 
are jet-dominate d (Falcke & Mai-koff 2003 iFalcke et al.1 120041; 
lYuan et al.l2002l : lFender et aL20 03). A compact, variable jet has 
in fact been detected in VLBI observations of M8 1 * (Bie tenholz 
et al., 2000; see also IVIarkoff et al., in preparation). iFalckd 
d 19961) demonstrate that a compact jet can explain the observed 
spectrum of M81* very well. In the jet model, the turnover fre- 
quency from a flat (or slightly inverted) radio spectrum to an 
optically thin power-law occurs at a frequency v,. In a simplified 
model this frequency depends on the jet power, Qj, and mass of 
the black hole, Mbh, as v, oc QJ^M^^^ dFalcke et alj|2004l) . In 



LLAGN the turnover occurs generally at (sub)mm wavelengths, 
while it is located in the near-infrared/optical regime for XRBs 
in the low/hard state (se e lCorbel & Fendeij|2002t iMarkoff et all 
120031; iFalckeet al.ll2004 . 

The simultaneous observations of M8 1 * at 3 and 1 mm pre- 
sented in this work are consistent with a turnover of the syn- 
chrotron emission from a jet in M81* into the optically thin 
part between 3 mm and 1 mm (with the exception of the first 
~4 hours of the 24 February light curve that would indicate a 
higher turnover f requency). This i s in g ood agreement with the 
data presented bv lReuter & Leschi ( Il996h that were, however, not 
acquired simultaneously. 
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Fig. 4. Spectral index of M81* on 24 February, 14/15 July, and 20 July 2005. The spectral index is defined hy Sy oc v", where Sv is 
the flux density at a given frequency v. The absolute uncertainties (see Tab. |2]l of the flux densities at 3 and 1 mm have been taken 
into account for this plot. Since we use upper limits on the absolute uncertainties, the plotted error bars are conservative. 



The evidence for the location of the turnover frequency be- 
tween 3 and 1 mm is ambiguous, however. The multi-frequency 
radio data from the coordinated campaign indicate that bumps 
are present in the radio spectrum (Markoff" et al., in prep.). These 
bumps are located at diff'erent frequencies at diff'erent epochs. 
Therefore, in spite of the presented mm-data, it may well be pos- 
sible that there is a submm-bump present in M8 1 *. Observations 
with the SMA at 345 Ghz that were obtained during the coordi- 
nated campaign on M81* indicate that the flux density of M81* 
increases toward the submm regime, in agreement with theo- 
retical predictions on the existence of a submm-bump (Markoff' 
et al., in preparation). Unfortunately, there was only one epoch 
(24 February 2005), where measurements with the PdBI and the 
SMA were simultaneous. 

For fine-tuning the models of the emission of M8 1 * it is im- 
portant to know the exact turnover frequency and whether and 
how it varies with activity of the source. Therefore there is an 
urgent need for more observational data for comparison with 
theory in order to understand sub-Eddington accretion and emis- 
sion. It is essential that the observations across the wavelength 
regimes are simultaneous. 

Due to the gap of several months there is probably no corre- 
lation between the February and July obse r vation s because the 
light curves at 2 cm presented by iHo et al.l (1 19991) indicate that 
radio outbursts of M81* hav e generally timescales < 5 months. 
Also, ISakamoto et all (|2001|) have found that M81* shows flux 
variations of factors < 2 at a wavelength of 3 mm on time scales 
< 10 days. However, the decaying light curve from July 19 may 
be related to the same event as the light curve on July 14. The 
spectral index during the two epochs is negative and of similar 
absolute value (see Tab. |3] and Fig. HJi. This behaviour is con- 
sistent with a decaying light curv e from a low-peaking flare in 
the model bv lValtaoia et al.l d 19921) that describes variability phe- 
nomena in AGN by a g rowth/decay of shocks in a jet (see also 
iMarscher&GeaJI 19851) . 

An estimate of the size scales of the relevant processes can 
be obtained from the time scales of the observed variability. The 
drop-off between 07:00 UT and 12:00 UT in the highest quality 
data set from 24 February 2006 shows that a > 5cr variability 
of the flux at 3 mm occurs on a time scal e of 5 hours (see als o 
the intraday variability at 3 mm found by Sakamoto et al.ll200Tl) . 
A similarly rapid change of the flux can be seen at 1 mm. Since 
no signal can travel faster than at the speed of light this cor- 
responds to an upper limit on the size of the source of merely 
~25 Schwarzschild radii, when a mass of 7.0 x 10^ is as- 
sumed for the black hole in M8 1 *. 



5. Summary 

We present three epochs of simultaneous 1 and 3 mm contin- 
uum observations of the LLAGN M81* that were obtained in 
the framework of a coordinated, multi-wavelength campaign 
(Markoff et al., in preparation). 

The observations of M8 1 * with the PdBI at mm-wavelengths 
confirm that M81* is a continuously varying radio contin- 
uum source as has been fou n d pre vi ously, mostly at larger 
wav elengths (see j HoetiD 1 1 9991; iBietenholz et al] 120001; 
Sakamot o et al.ll200Tl) . The measurements present the first un- 
ambiguous detection of M8 1 * at 1 mm and moreover show that 
the source is continuously variable at this wavelength as well. 
The amplitude of the variability is observed to be generally 
larger at 1mm than at 3 mm by a factor of rough ly 1.5. This 
agrees well with the trend found bv lHo et aTl ( Il999l) that the am- 
plitude of the variability increases with frequency in LLAGN. A 
similar behaviour is found in the far weaker source Sgr A* (e.g. 
iHerrnstein et al.ll20 04'; Miva zaki et"aLll2004 . 

The shortest variability time scales of our observations 
give upper limits on the size of the emitting region of 
~25 Schwarzschild radii, assuming a black hole mass of 7.0 x 
10^ Mq. 

The decaying light curves observed on 14/15 and on 20 July 
2005 may have been related to the same radio out burst. They are 
consistent with a generalized shock-in jet model (IValtaoia et al.l 
.19921) . 

The simultaneous measurements of M8 1 * at 3 and 1 mm 
are consistent with a turnover of the flux in the mm-to-submm 
regime as predicted by models for a jet dominated source. Some 
ambiguity remains concerning the exact peak frequency and the 
related interpretation of variability events. 

The observations confirm previous findings that there are 
many similarities between M81* and Sgr A*, the source re- 
lated to the supermassive black hole at the center of the Milky 
Way. This underlines the similarity between LLAGN despite of 
several orders of magnitudes of difference between their lumi- 
nosities. M81* can serve as a bridge from the extremely sub- 
Eddington Sgr A* toward higher luminosity LLAGN. The radio 
emission from M81* is apparently dominated by a compact jet. 
The same may be the case for Sgr A*, although such a jet has 
not yet been discovered due to the strong interstellar scattering 
toward this source. This work poses a clear case for follow-up 
observations of the highly sub-Eddington nucleus M8 1 * and for 
the necessity of simultaneous multi-wavelength observations in 
general in order to constrain the physical mechanisms in black 
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holes that accrete at rates several orders of magnitude below the 
Eddington limit. 
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